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in the upper right, with the practical absence of analyses in the upper 
left and lower right quadrants, is very evident. 

The chemical composition of igneous rocks is so complex and the 
number of their chemical constituents is so great; the conditions of 
differentiation and solidification are so varied and so complicated; and 
the knowledge of the application of the laws of physical chemistry to 
them is as yet so meager by reason of the lack of data, that some irregu- 
larities and apparent exceptions are to be expected. Furthermore there 
are other correlative relations which have not been dealt with here, and 
it is conceivable or indeed probable that one or more of these may under 
certain conditions, whether of composition, differentiation or solidifica- 
tion, supersede the one which is the subject of the present discussion 
without invalidating the truth of this. 

It would be premature to enter here into a discussion of the causes 
of the correlation of these four elements. Any such would necessitate 
the consideration of the other correlations which have been detected. 
It may be said, however, that even though such concomitant variations 
may be due to similar solubility relations, or some such physico- 
chemical factor, yet that such similarities are themselves presumably 
due to certain intimate relationships between the elements which are 
commonly and somewhat loosely called affinities. 

1 J. H. L. Vogt., 2s. prdkt. Geol., 1898, 326. 

2 J. F. Kemp., Ore Deposits, 34-37 (1900). 

3 L. de Launay, La Science Gtologique, 637, 190S; Giles Mineraux, 1, 46-51 (1913). 
* W. F. Hillebrand., Bull. U. S. Geol. Sun., 305, 21 (1907). 

6 H. S. Washington., Trans. Amer. Inst. Min. Eng., 1908, 749-767. 
' H. S. Washington., Prof. Paper U. S. Geol. Sun., 14, 1903. 
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In earlier papers I have considered the effect of a linear transfor- 
mation of dependent variables upon the solutions of ordinary linear dif- 
ferential equations in the vicinity of a singular point. 

This type of transformation led me to the notion of equivalence which 
is fundamentally important for the classification of singular points. 

Ordinary linear differential equations also preserve their form under 
an arbitrary transformation of the independent variable. I shall prove 
here that this second type has no additional significance for the purposes 
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of classification since the solutions are carried over into the solutions 
of equivalent differential equations by such transformations. 1 The 
precise statement of the result is contained in the following 

Theorem. Let Y (*) be a matrix 2 of functions y tj (x) (i, j = 1, 2, • • • , n) , 
each analytic for \x\^.r, with determinant Y (x) not zero for \x\^r. 
A necessary and sufficient condition that Y (x) forms the matrix solution 
of a linear differential system 

d JL = A{x)Y, 
dx 

where the elements a i} - (x) of A (x) are analytic (q = 0) or have a pole of at 
most the qth order (q>0) at x = <&, is that, for every function <p (x) 
of the form 

<p (x) = x + I x~ q + m x ~ 9-1 + • • • , 

the matrix M (x) defined by the matrix equation 

Y 0)) =M(x) Y(x) 
is composed of elements m tj (x) analytic at x= <». 

Let us prove first that the condition is necessary. 

The function M(x) defined for any choice of ip(x) is single-valued as 
well as analytic for \x\ sufficiently large. In fact, when x makes a cir- 
cuit of x= co , Y(x) is altered to Y(x)C where C is a matrix of constants 
of determinant not zero. But since ip(x) —x remains finite as x becomes 
infinite, <p (x) will make essentially the same circuit if \x\ be sufficiently 
great. Thus Y((p(x) ) changes to Y(<p(x) )C, and Mix) is unaltered in 
value. 

Now the expression for M (x) 

Y Oq) F _1 OJ (x = <p(x);x m = x) 

may also be written as a matrix product 

[Y (*„) F- 1 (*,)] [Y ( Xl ) F- 1 (x 2 )]...[Y (*„_,) F- 1 (*„)]. 

Here we will assume that the values x , x h x%, . . . , x m of x are 
equally spaced, each after the first differing from its predecessor by 
Ax = (xi—Xm)/m. 

From the differential equation itself, it is clear that 

Y (Xi) - Y (x i+1 ) = Ax[A (x t ) + E { ] Y (x i+1 ) (i = zero, 1, • • • , m - 1) 

where Ei will have arbitrarily small elements if Ax is small enough. 
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For |*| taken large enough, it is also clear that we will have an equality 
of the form 

| a (x) | < a \x q | (i,j, = 1, 2, • • ■ , n) 

for the elements of A (x) . Hence if / denotes the matrix with elements 
8ij(Sij = 0, »=fj-; S»» = l), and if A denotes the matrix whose elements are 
all equal to the positive constant a, we find that the elements of the 
matrix 

Y (*,-) Y- 1 (x i+1 ) or / + Ax (A fe) + E t ) 

are less in absolute value than the corresponding elements of the matrix 

I + 2\&x\A\xi\". 
Also on account of the form of <p(x) we have for |*| large enough 
| x — x m | g k | x t I ~ q , (k independent of i or x). 
and thence 

|Ax|S — \xi\-". 
m 

Substituting this value of A* above, we conclude that the elements 
of Y(xi) Y'^Xi+i) are less in absolute value than the corresponding 
elements of 

I+*A. 

m 

Thus, by the product formula, Mix) will have elements whose absolute 
value is less than that of the corresponding elements of 



(i + —a\ 

\ ml 



whose elements are obviously limited in absolute value for all values of m. 

Consequently the elements of M (x) are limited in absolute value for 
|*| large enough, and are necessarily analytic at infinity. 

Let us show now that the stated condition is sufficient. 

To this end we will consider the family of functions 

<p (x) =x + px~ g 

where p is a parameter. 

Defining M(x,p) by the equation 

Y (x + P x~ q ) =M(x, p) Y(x), 
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we observe that M (x,p) is made up of elements single-valued and ana- 
lytic in x and p if 

|*|^ 1 +r, |p|^l, 
inasmuch as in this event 

Itf + pz-'I^M -\x~ g \^\x\-l^r 

Moreover, in view of the stated condition, the elements triij(x,p) 
of M(x,p) are analytic in x at x = oo . 
Consider now the integral formula 

I^fcp)--^ fa^ (|.|>l + r) 

x 2w v — 1 ^ £ (£ — *) 

where the integration is performed around the fixed circle C: |f| = l+r 
in a positive sense. This reduces to the Cauchy integral formula by 
the substitution x=l/x', £=l/£', which holds since m(x,p) is analytic 
for x ^ l+r and at x= oo. 

We now perceive from the form of the integrand that ma (x,p) must 
be analytic in x and p for x = oo . 

Differentiating the equation of definition for M(x,p) with respect 
to p we obtain 

x-*4- Y(x + px-') = 4-M (x, p) Y(x). 
ax dp 

The elements of dM(x,p)/dp are also analytic in x and p at *= oo . Let- 
ting p approach zero we obtain, for \x\ > l+r, 

x -" d I^l = ±M(x,0)Y(x) 
dx dp 

which establishes at once that Y(x) is the solution of a differential 
system of the desired form. 

1 This result was partly known to me in 1908. For a special case, see Trans. Amer. Math. 
Soc, 14 462-476 (1913), in particular pp. 475-476. 

2 For the elements of the theory of matrices here used see Schlesinger, VorUsungen ilber 
linearen Differentialgleichungen, pp. 18-19. 



